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Abstract 

We recently introduced a wave-mechanical model for quasi-elastic neutron scattering (QENS) 
in proteins. We call the model ELM for “Energy Landscape Model”. We postulate that the 
spectrum of the scattered neutrons consists of lines of natural width shifted from the center by 
fluctuations. ELM is based on two facts; Neutrons are wave packets; proteins have low-lying 
substates that form the free-energy landscape (EEL). Experiments suggest that the wave packets 
are a few hundred micrometers long. The interaction between the neutron and a proton in the 
protein takes place during the transit of the wave packet. The wave packet exerts the force F{t) = 
dQ{t)/dt on the protein moiety, a part of the protein surrounding the struck proton. Q{t) is 
the wave vector (momentum) transferred by the neutron wave packet to the proton during the 
transit. The ensuing energy is stored in the energy landscape and returned to the neutron as the 
wave packet exits. Kinetic energy thus is changed into potential energy and back. The interaction 
energy is proportional to Q, not to Q^. To develop and check the ELM, we use published work on 
dehydrated proteins after reversing improper normalizations. In such proteins only vibrations are 
active and the effects caused by the neutron momentum can be studied undisturbed by external 
fluctuations. ELM has predictive power. Eor example it quantitatively predicts the observed 
inelastic incoherent fraction S{Q,T) over a broad range of temperature and momentum Q with 
one coefficient if 5(0, T) is known. 
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I. BACKGROUND 


The goal of our research is the exploration of protein dynamics. Life can only exist if 
its building blocks can move [T]. A vast arsenal of techniques has been brought to bear on 
protein dynamics. Our own work started many years ago using the Mossbauer effect and 
flash photolysis. Recently we began to look at neutron scattering [2]. In these experiments 
neutrons with wave vectors q hit the protein. The intensities S{Q,T) of the scattered 
neutrons are measured as a function of the temperature T at different scattering angles, 
characterized by their wave vectors q'. The wave vectors q and q' determine the transferred 
wave vector, Q = q ~ q'■ We use the standard notation where the wave vector q is related 
to the momentum by p = hq. For convenience we will call Q also “momentum transfer” or 
simply momentum. Neutron scattering is treated in excellent texts and reviews [3H7]. The 
salient experimental observation is easy to describe: The scattered slow neutrons appear to 
show an elastic line and a broad quasi-elastic band. The currently accepted model explains 
the quasi-elastic band as being due to the Doppler effect caused by spatial motions of the 
target; we call it SMM for “spatial motion model”. In the SMM the spectrum of the scattered 
neutrons is separated into two components, an elastic line and a quasi-elastic band. The 
quasi-elastic band is assumed to be composed of broad Lorentzians centered at the elastic 
line. In the simplest model, the quasielastic Lorentzians are attributed to motions of the 
target proton in a Gaussian well with mean-square displacement (MSD) (x^). S{Q,T) is 
then given by the Debye-Waller (Lamb-Mossbauer) relation 

S(Q,r)=exp{-Q={i"(r)>} (1) 

Here (x^(T)) is the average MSD of the protons in the protein. If Eq. does not fit the 
experimental data, “non-Gaussianity” and “dynamic heterogeneity” [HHn] are invoked. 

A number of issues caused us to look critically at this model: We could not hnd a 
transparent description of the quasi-elastic neutron scattering (QENS) and of the physics 
behind the momentum transfer Q. The separation of the spectrum into an elastic line and 
a quasielastic band made no sense to us. The existence of broad homogeneous Lorentzians 
centered at the elastic line is difficult to accept because the uncertainty relation would im¬ 
ply very short lifetimes. Eq. [^predetermines the explanation of the quasi-elastic scattering: 
QENS must be due to motions in real space to account for the large mean-square displace¬ 
ment (x^). The focus on spatial motions may have been influenced by the early scattering 
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theories that assumed the absence of low-lying substates [12]. Non-Gaussianity and dynam¬ 
ical heterogeneity appear like hrst-aid patches to save the model. These issues caused us to 
abandon the Gaussian approximation and its generalizations and search for a physics-based 
model with a testable role for Q. The hrst step in this endeavor led to a model that does 
not involve Eq. [Tjand that is based on the free-energy landscape concept. We call it “energy 
landscape model” or ELM HSllll]. It answers some questions, as we describe in the next 
section. The meaning of Q, however, remained unclear. We propose to explain Q here. We 
therefore consider this paper part II of the PNAS publication |14j . 

II. THE ENERGY LANDSCAPE MODEL (ELM) 

The energy landscape model is based on two pillars: Neutrons are de Broglie wave 
packets [T51 ITB] : proteins have a large number of low-lying conformational substates (GS) 
that are organized in the free-energy landscape (EEL) [13 [IE] • The wave nature of the neu¬ 
tron implies that the neutron-proton interaction time can be of the order of hundreds 
of ns. The existence of the GS implies that proteins can change their conformation and 
their energy during the passage of the wave packet. The two pillars permit the construc¬ 
tion of the ELM as we have described previously [131 IH|. Fig. describes the essential 
properties of the EEL. Fig. [^A) is a schematic of the density D(E) of the substates in 
the EEL. D{E)dE is the fraction of substates with energies between E and E -|- dE. We 
assume that D{E) is constant below about 50 meV. The population p{E) of the substates 
depends of course on the temperature. In Fig. [^B) a few GS have been selected. Gonsider 
a protein in a particular substate. In the absence of fluctuations it stays in this substate for 
a long time and neutron scattering yields a single line at AEq = 0. The width of the line 
is given by the resolution function of the instrument. Proteins have a broad spectrum of 
thermal fluctuations. Below about lOK, the fluctuations are too slow to induce observable 
transitions among substates. As the temperature is increased, vibrations become faster and 
their density increases [T^ [20] . The vibrations switch a protein from GS to GS, with average 
fluctuation time Tf. Assume that the system is observed for a time r„. If is much shorter 
than Tf no temporal information can be recovered. However, if is much longer than Tf the 
energy shift during the time Tn can be extracted as is suggested by Fig. [^B). The protein 
fluctuates approximately r„/r/ times during the observation time r„, each time jumping 
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from substate to substate and winding up in a substate shifted by the energy IS.E. The 
intensity S'(0,T) of the peak decreases and a Lorentzian appears aX E + AE. The shifts 
AE can be positive or negative, and the result is an inhomogeneous spectrum, composed of 
shifted Lorentzians as sketched in Fig. [^C) The spectrum is dynamic; there is no difference 
between the central line and lines in the quasi-elastic band. No substate is permanently 
occupied; any protein can transiently be in any substate. There is no separation into an 
elastic line and a quasielastic band. All individual spectral lines have the same width, of the 
order of a few neV. The energy of the initial substate is irrelevant for AE, only the change 
in energy during the time counts. Different proteins experience different AE and the en¬ 
semble produces a broad spectrum as in Fig. [^C). In neutron scattering the relevant time 
Tn is determined by the optics in the neutron source na and can be of the order of hundreds 
of ns. In the Mossbauer effect, the lifetime tmo =141 ns replaces r„. In the Mossbauer effect 
the natural line width is extremely small and the decrease of the central intensity and the 
resulting line broadening can be observed already below 20 K [21] . In neutron scattering, the 
observed line broadening is influenced strongly by the resolution function of the instrument 
and occurs at higher temperature than in the Mossbauer effect. 

Vibrations are not the only fluctuations that produce transitions among substates. In 
hydrated proteins the a-fluctnations in the solvent and the /3h-fluctnations in the hydration 
shell can cause protein motions [221126] . QENS studies of systems involving more than one 
type of fluctuations are harder to interpret and hide the simplicity of the ELM. We therefore 
use in the present paper only dehydrated proteins, where vibrations rule. 


III. FORWARD INCOHERENT NEUTRON SCATTERING 

In principle, the fraction S'(0,T) can provide crucial information. In practice, however, it 
is difficult to measure it directly. S'(0,T) must thus be determined indirectly. One approach 
is to average a number of measurements with small Q values as done by Hong et al.|27]. 
This approach yields reasonable data, but for Q > 0. We have used the data at different Q 
values in the supplementary information to the paper by Hong et al. to extrapolate S{Q, T) 
to Q = 0. The result, after normalizing S'(0,0) to one, is given in Fig. Also shown in 
Fig. I^is S'(0.8A T) to indicate that there is a measurable difference between S'(0.8A T) 
and 5(0, T). The reason why S{Q > 0,T) is smaller than 5(0, T) will become clear later. 
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Fig. 1^ also suggests that S{0,T) decreases little below about 100 K. This effect is due to 
the broad resolution function in neutron scattering. The decrease in 5(0, T) becomes only 
observable when a sizable fraction of proteins have moved away from the center. However, 
comparison with Mossbauer data izn, where the energy resolution is much better, indicates 
that the decrease below 100 K is indeed small so that Fig. [^is a good approximation. 

IV. DENORMALIZATION 

Normalization is the migraine headache of QENS. Improper normalization can warp the 
data. The correct normalization is 5(0, 0) = 1. We did not find data approaching this ideal. 
Most data are “normalized” using one or sometimes two procedures: 

i. The assumption that all motions are absent at very low temperatures, for instance at 
To = 10K, leads to the normalization 5(Q,T)/5(Q,To). This procedure is incorrect 
because quasi-elastic effects exist even at 10 K as we will show. Denormalization is 
done by multiplying the published data by S{Q,Tq) to recover S{Q,T). 

ii. Assuming that the elastic effects can be removed without changing the quasi-elastic 
shape leads to division by the elastic fraction, 5(Q,T)/5(0,T). This normalization 
changes the shape of the published S{Q, T). Denormalization is in principle performed 
by multiplying the published 5(Q,T) by the elastic fraction 5(0,T). If the 5(0,T) 
used in normalization is not given in the publications, we denormalize by using 5(0, T) 
from Fig. While this procedure may lead to some deviation from the true value, it 
corrects the data in the right direction. 

As an example of denormalization we select QENS data by Wanderlingh and collaborators 
[28] . These authors measured incoherent quasielastic spectra of glutathione S-transferase 
(GST). Their spectrum, after digitizing and smoothing, is given in Fig. [^A). Denormal¬ 
ization i is not needed, because the authors used a vanadium standard to get the correct 
S{Q,T) and did not divide by S{Q,T < 30K). Denormalization ii is required because the 
data were normalized to S{Q,T) = 1. Because 5(0, T) is not given in their publication we 
use 5(0, T) from Fig. The denormalized S{Q, T) is plotted as function of Q in Fig. [^B). 
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V. THE MOMENTUM TRANSFER Q EXCITES THE PROTON MOIETY 


The data of Wanderlingh et ah in Fig. |^B) expose the surprising fact that the fraction 
S{Q,T) depends strongly on Q even at 30 K. For instance ^(dA \ 30iF) is only 0.76. 
Motions at 30 K are only due to thermal vibrations. These decrease S'(0,T) by at most a 
few percent at 30 K as shown in Fig. Fig. [^also shows that S'(0, T) = 0.76 corresponds to 
a temperature T* ^ 400K. Thus the momentum transfer Q = dA ^ increases the effective 
temperature by AT*(dA ^ 370K). Calculation of AT*(Q) for a number of momenta leads 

to Fig. I^C). This figure shows that AT*(Q) is proportional to Q and suggests the relation 

ataq) = xg (2) 

Here y is a coehcient with a value of ~ 88 K A ^. 

T* cannot be a real temperature. A heat spike would dissipate in ps. A simple model 
suggests that the energy AE* k. ksAT* can be stored as potential energy in the FEL. 
Consider an ideal elastic spring with spring coefficient k. If a force F compresses the spring 
by the distance x the spring reversibly stores the energy In proteins, the basic features 

may be similar. The force could be due to the neutron wave packet and the moiety around 
the proton could act like an elastic spring. While a change in temperature would be shared 
by the entire sample, a structural change can be local. It could for instance correspond to a 
small change in position of a side chain close to the target proton. The excited moiety thus 
could be the mechanical analog of a thermal hotspot. 

The suggested mechanism for the quasi-elastic neutron scattering is drawn in Fig. In 
Fig. I^A) the wave packet, approximated by a Gaussian, hits the proton with the momentum 
Q{t). The impact creates a force F{t) given by. 

Fit) = dQ{t)/dt. (3) 

Fit) is sketched in Fig. |^B). During the hrst phase of the transit. Fit) is positive and the 
wave packet compresses the moiety by a distance x. This compression increases the potential 
energy AE*it). During the second phase Fit) is negative and the energy is returned to the 
wave packet. The protein exits with essentially the full energy, modulated by the quasielastic 
fluctuations that produce the broad band. The stored energy is sketched in Fig. |^C). The 
form of AE*it) depends on the compressibility of the moiety. In the simplest case the process 
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leads to AE* ^ x^. In proteins the compressibility decreases markedly upon compression 
[2n] and the simplest relation may not apply. Fig. |^C) thus only gives a guess as to the 
form of the stored energy AE*{t) as a function of the transit time. 

The energy input due to the momentum transfer could also be stored as kinetic energy 
in local vibrations. 

Support for an elastic spring model comes from the fact that polymers can store energy 
|30j . Shock wave experiments on polymer solids also involve energy storage [3TH35] . In an 
experiment by Dlott and coworkers a 50 pm thick aluminum plate hits a solid polymer target 
with a velocity of 1 km/s. The shock compresses the polymer solid. The deposited energy 
excites the FEL of the polymer, leading to an increase in the local temperature by about 
100 K [36]. 

The observation that AT{Q) in neutron scattering is of the same order as seen in the shock 
wave experiments permits an estimate of the size of the target in the neutron scattering. 
In the shockwave experiment the momentum density is about Pg = 10^kgm“^s“^. In 
the QENS experiment the relevant momentum is about 10“^^kgms“^. We 

interpret the ratio Pn/Ps ~ 10“^^ m^ as the cross section for the neutron-proton scattering. 
It is remarkably close to the experimental value |37| of 2.4 x 10“^^ m^ and supports the 
assumption that the neutron interacts with a single proton. 

The comparison between the shock wave and the QENS processes cannot be taken too far. 
In the shock wave experiment, the aluminum plate is destroyed and the energy is ultimately 
dissipated as heat. In neutron scattering the wave packet interacts with a single proton 
that then puts pressure on its immediate surrounding. The pressure in turn changes the 
conformation of the proton moiety [38] thereby storing the energy reversibly in the FEL. 
The interaction is local and the FEL is only changed locally. The energy is not lost but 
returned to the neutron. 

VI. THE EFFECT OF VIBRATIONS 

Conformational motions can be caused by thermal vibrations, /3/i-fluctuations in the hy¬ 
dration shell, or a-fluctuations in the environment. If they are simultaneously present it 
is difficult to unravel their effects. We consequently selected a dehydrated protein where 
thermal vibrations dominate. We use QENS data of Nakagawa, Kamikubo, and Kataoka 
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on dehydrated Staphylococcal nuclease (SNase) [39]. Fig. [^A) presents some of their data 
after digitizing and smoothing. The data imply that both normalizations i and ii have 
been used. Denormalization i is straightforward; the data are multiplied by S{^,T) from 
Fig-i Denormalization ii is more involved because the measured S{Q, lOK) has not been 
published. However, the data in Fig. [^B), together with Fig. permit one to reconstruct 
S{Q^ lOK). We assume that Eq. holds even at the lowest temperature. Fig. [^C) predicts 
that the temperature increase for instance at Q = dA ^ is about 352 K and T* = 362K. Fig. 
[^then gives S{4:A lOK) ^ 0.79. Thus values for S{Q, T) can be found even for T = lOK, 
and denormalization ii can be performed. Fig. j^B) gives the denormalized S{Q,T). 

Information on protein dynamics is contained in the temperature dependence of S{Q, T). 
This dependence is difficult to see in Fig. |^B). We therefore redraw S{Q,T) as function of 
temperature in Fig. [^C). Also shown in Fig. 5(C) is S{0,T), taken from Fig. The result 
is surprisingly simple: The S{Q,T) curves have essentially the same shape as the elastic 
S{0,T), but are shifted to lower temperatures. This shift is at hrst counterintuitive. Why 
should the energy input due to the neutron lead to lower temperatures? The answer is in the 
difference between T and T*. T is the ambient temperature of the entire sample. T* gives 
the energy of the excited state in the FEL produced by the pressure of the neutron wave 
packet. T* is higher than T by AT*{Q) given in Eq. Thus the protein moiety measured 
at the low temperature T is actually at the temperature (energy) T*. 

To quantify the shifts between T and T* we plot in Fig. |^D) the temperatures of S{Q, T) 
at hxed values of S{Q,T). The curves are parallel and show that for a hxed S(Q, T) the 
shift is proportional to Q and given by Eq. |^with y = 88KA. The neutron increases the 
temperature (energy) of the proton moiety by AT*{Q) = xQ, or 

A(Q,n = A(0,T + yg). (4) 

The incoherent elastic fraction S{0,T) and the coefficient y together predict S{Q,T) for 
all values of Q. S{0,T) with dehydrated proteins is given by the spectrum of thermal 
fluctuations. The only £t parameter is the coefficient y in Eq. 


VII. CONCLUSIONS 


Incoherent quasi-elastic neutron scattering (QENS) has traditionally been described as 
follows: Above a certain temperature, the sharp elastic line decreases in intensity and a broad 
“quasi-elastic” band, centered at the elastic line, appears. In the conventional approach 
to neutron scattering, used for more than fifty years, the two phenomena are ascribed to 
spatial motions, for instance spatial diffusion M. No clear interpretation of the effect of 
the momentum transfer Q is given. 

In ELM, the energy landscape model, the central line and the broad band are part of 
one spectrum, consisting of lines with essentially the natural width, and shifted by fluctua¬ 
tions. The neutron is a wave packet, a few hundred pm long. During the transit, the wave 
packet exerts the force F{t) = dQ{t)/dt on the proton moiety. The ensuing energy AE*{t) 
is temporarily stored in the free energy landscape. Thus the kinetic energy of the neutron is 
transformed into potential energy, not heat. Some or nearly all of the energy is returned to 
the outgoing wave packet. The scattered neutron has essentially the full energy. During the 
passage vibrations induce transitions among the substates of the EEL, thereby producing 
the quasielastic band. To test the model, we have used published data after denormaliza¬ 
tion. Despite the fact that ELM explains the data, the currently accepted model must also 
contribute. Fluctuations in the energy landscape are connected to conformational motions 
in real space. A future theory must connect the two views. 

The ELM solves or avoids the problems that we listed in Section The ELM provides a 
physics-based explanation of QENS. The central line and the broad band are components 
of the same spectrum. The components are narrow and do not cause problems with the 
uncertainty relation. Non-Gaussianity and dynamic heterogeneity are not needed. The role 
of Q is clear; it creates the force that the passing neutron wave packet exerts on the proton 
and that leads to the excitation of the protein moiety. The effect of the wave packet is 
proportional to Q, not to Q^. Spatial diffusion, which is difficult to explain at very low 
temperatures, does not enter. The observed spectra with non-zero Q are quantitatively 
predicted once S'(0,T) and the coefficient y are known. The simplicity of the model is 
bought, however, with new problems. Can the neutron wave packet slow down and then 
speed up in its passage? Can the force due to the neutron wave packet reversibly excite the 
proton moiety without dissipation? What is the nature of the coefficient y? 
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Finally a remark about the theoretical underpinning of the interaction of the wave packet 
with the proton. How can the long-lasting interaction be described by the instantaneous 
hrst Born-approximation? Baym and Pethick have treated this problem by assuming an in¬ 
stantaneous contact interaction and replacing the customary plane waves with wave packets 

m- 
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FIG. 1. (A) Hypothetical density of substates D(E). (B) Random walk of a protein in the energy 
landscape, caused by thermal fluctuations. Tf is a typical vibration frequency, 6E typical jump 
energy between substates, and AE is the final substate energy shift. (C) AE has shifted a line 
from it initial position in the center line to its hnal position. 
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FIG. 2. The elastic fractions S{Q,T) for Q = OA ^ and Q = O.sA ^ for dry green fluorescent 
protein (GFP). Data taken from Figure SI of the Supplementary Material of Hong et ah |27j . The 
curve for Q = 0 was obtained by extrapolating the data for three wave vector transfers to Q = 0. 
The dashed line results from an extrapolation between 250 and 300 K to 400 K. 
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FIG. 3. (A) log{S{Q,T)) vs. for glutathione transferase (GST) at hydration of 0.5 and tem¬ 
peratures 30, 105, and 165 K after smoothing of digitized data. Adapted from Wanderlingh et 
al. [28] where the data is normalized to 5(0, T) = 1.0. (B) S{Q,T) plotted versus Q for the 
de-normalized data using Fig. See text. (C) The relation AT{Q) = xQ from Fig. [^and[^-B 
as described in text, x = 88K A. Smoothing based on Lowess (Weighted) algorithm implemented 
in KaleidaGraph 4.5, Synergy Software 
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FIG. 4. (A) Q{t) is the intensity of the neutron wave packet at the proton at the time t. (B) The 
force F{t) exerted by the wave packet on the proton at time t. (C) The energy AE*{t) deposited 
in the proton moiety at the time t. 
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FIG. 5. (A) Normalized S{Q,T) vs Q for dry staphylococcal nuclease (SNase). Dashed lines 
show linear extrapolation from smallest measured Q to Q = 0. Adapted from Nakagawa et al 
[39]. (B) S{Q,T) vs Q for dry SNase after de-normalizations i and ii. Dashed lines are linear 
extrapolations from smallest measured Q to Q = 0. Simulations of S{Q,T) at T = 0.1 and 10 K 
(all dashed/dotted ) as explained in text. (C) S{Q,T) vs T for dry SNase including the curve for 
Q = 0 (dotted). Smoothing in (A), (B), and (C) is based on the Stineman (Smooth) and Lowess 
(Weighted) algorithms implemented in KaleidaGraph 4.5, Synergy Software, 2013. (D) Fig. 
shows that the intensities S{Q,T) at different values of Q have the same overall temperature 
dependence, but are shifted to lower temperatures with increasing Q. To quantify the shifts, Fig. 

presents cross sections through the data in Fig. at selected values of S{Q,T) from 0.3 to 
0.9. The cross sections are parallel with y ~ 88KA as for glutathione S-transferase. 
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